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LINEAR  AND  NONLINEAR  PULSE  PROPAGATION  IN  OPTICAL  WAVEGUIDES 
1.  INTRODUCTION 

The  possibility  of  transmitting  undistorted  pulses  with  high  peak 
powers  guided  by  fibers  could  prove  useful  and  find  application  in 
various  fields.  Intense,  ultrashort  light  pulses  guided  by  fibers  could 
be  used  to  achieve  ultrahigh  data  rates  in  communications,  to  perform 
surgery  in  medicine,  for  cutting  and  welding  in  industrial  processes, 
etc.  The  nonlinearity  which  governs  the  propagation  characteristics 
arises  from  the  real  part  of  the  guide  refractive  index,  although 
dissipative  nonlinear  processes  could  in  principle  affect  the  pulse  as 
well. 

The  theoretical  problem  of  transmission  of  nonlinear  pulses  in 
optical  waveguides  is  both  interesting  and  challenging,  since  it  involves 
the  simultaneous  interplay  of  transverse  confinement,  dispersion  and 
nonlinearity.  Solutions  in  many  areas  of  physics  have  already  been 
obtained  (SCOTT,  A. C.  et  al.,  1973,  for  example)  and  considerable 
knowledge  has  been  gained  for  the  one-dimensional  pulse  propagation  case, 
while  two-  or  three-dimensional  nonlinear  propagation  is  still  in  its 
infancy  stage  (GERSTEN, J. I. ,  et  al.,  1975?  KAW,  P.K.,  et  al.,  1975? 
WILCOX,  J.Z.,  et  al.,  1975?  ZAKHROV, ,  V.E.,  et  al.,  1979). 

A  critical  limitation  in  the  realization  of  the  full-bandwidth 
capability  of  optical  transmission  systems  is  pulse  broadening  resulting 
from  dispersion  effects  (GLOGE,  C. ,  1971?  KECK,  D.,  1976,  for  example). 
Although  dispersion  can  be  minimized  by  an  appropriate  choice  of  guide 
materials,  geometry  and  operating  frequency,  it  is  apparent  that 
dispersion  effects  will  be  detrimental  for  sufficiently  narrow 
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pulsewidths.  To  overcome  this  problem  the  nonlinear  dependence  of  the 
refractive  index  on  pulse  intensity  may  be  used.  Here  the  dispersion 
effects  result  in  the  broadening  of  the  pulses,  while  nonlinearity  tends 
to  sharpen  it.  It  is  the  appropriate  interplay  of  these  opposite  effects 
which  can  lead  to  a  stable  soliton  solution  for  the  optical  pulse.  But 
nonlinearity  can  also  introduce  a  host  of  other  effects  (CHIAO,  R. Y. ,  et 
al.,  1964;  HASEGAWA,  A.,  et  al.,  1973;  TZOAR,  N. ,  et  al . ,  1976;  SHIMIZU, 
F. ,  1967)  such  as  self -focusing  and  self-phase  modulation.  However,  each 
of  these  effects  is  dominant  only  in  certain  regimes  characterized  by  the 
dispersion,  peak  intensity  and  nonlinear  parameters. 

Self -focusing ,  which  occurs  for  finite-sized  beams,  requires  a  very 
high  critical  intensity.  This  phenomenon  is  an  instability  which  causes 
a  finite-sized  beam  to  collapse  to  a  point,  due  to  an  effective  focusing 
lens  induced  by  its  own  intensity.  The  condition  for  self -focus  is 
independent  of  the  dispersion  parameter  and  represents  a  competition 
between  nonlinearity  and  diffraction  effects.  Indeed,  for  Gaussian  beams 
the  condition  for  self -focusing  to  occur  is  given  by 

(II 

Here  k=  2^^ /^,  where^  is  the  free  space  wavelength:  is  the 

linear  (nonlinear)  dielectric  function;  and  and  F^are,  respectively, 
the  radius  and  on-axis  field  intensity  of  the  beam.  This  condition  is 
not  affected  by  the  guide’s  parameters.  The  guide,  by  balancing 
diffraction,  provides  stable  trapping  of  the  beam  at  a  finite  size  even 
for  fields  below  the  critical  field  (the  critical  field  is  defined  by  Eq. 
(1)  when  the  equal  sign  is  taken) .  For  glasses  one  finds  the  critical 
power  for  self  focusing  to  be  of  the  order  of  10^-  10®  W/cm^,  much 
greater  than  the  power  needed  for  soliton  propagation. 
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The  soliton  problem  in  optical  waveguides  appears  to  have  been 
considered  first  by  HASEGAWA  and  TAPPERT  (HT)  (1973)  but  these  authors 
did  not  take  detailed  account  of  the  transverse  inhomogeneity  in  their 
analysis.  Subsequently,  JAIN  and  TZOAR  (JT)  (1978)  introduced  an 
approach  that  accounted  for  the  transverse  inhomogeneity  of  the  waveguide 
in  an  average  fashion,  and  were  thereby  able  to  demonstrate  that  in  fact 
typical  fiberguides  could  indeed  support  "bright"  solitons.  This 
contrasts  with  the  results  of  HT  who  found  that  anomalous  material 
dispersion  or  very  large  core-clad  differences  (^n'^0.5)  are  required 
for  "bright"  solitons  to  exist.  The  approach  of  JT  has  been  generalized 
by  BENDOW,  GIANINO,  TZOAR  and  JAIN  (BGTJ)  (1980)  to  include  variational 
analysis  and  to  treat  longitudinal  inhomogeneities,  and  by  CROSIGNANI, 
PAPAS  AND  DI  PORTO  (1980)  to  consider  the  role  of  intensity  fluctuations 
on  nonlinear  soliton  propagation  in  optical  waveguides. 

2.  THEORY  OF  PULSE  PROPAGATION 


Consider  an  optical  waveguide  having  a  dispersive  and  nonlinear 
refractive  index  of  the  form  _ 

e)  -  nt,  (n,M)  + 


-12. 


(2) 


where W  is  the  frequency,  E  is  the  electric  field  and  n2  the  nonlinear 
coefficient,  taken  to  be  independent  of  frequency.  Moreover,  n^  is 


represented  by  its  local  approximation,  i.e. 

(ri,U))  ■=  ^  (it) 


(3) 


In  what  follows  we  consider  propagation  of  pulses  which  are  narrowly 
centered  about  a  given  frequency  and  assume  that  nj^(W)  is  a  slowly 
varying  function  of  (a)  in  the  vicinity  of  .  The  electric  field  is 
taken  to  be  .  :  (at  ^  /•>  4-) 


E  (i^A) 


(4) 
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Where  e  is  a  unit  vector  in  the  transverse  direction,  r=  (^,z),  and  z  is 
the  direction  of  propagation.  We  then  obtain  for  A  the  wave  equation: 


[Vf*  ^  {(Z)  l<'.  f(ii) 

-  ( fe,  kl  +  ^  I  Ar/C  =  0 


(5) 


In  Eq.  (5)  k(tJ)  =nj^((»i)W/c  and  the  primes  indicate  derivatives  with 
respect  to  CJ  at  U  =U)^.  The  appearance  of  the  time  and  frequency 
derivative  only  to  second  order  is  a  consequence  of  using  the  slowly 
varying  envelope  approximation  in  deriving  Eq.(5). 

For  the  waveguide  problem  f  is  always  inhomogeneous  in  the  transverse 
direction,  i.e.,  f(r)  =  f  (f).  Under  these  conditions  it  is  extremely 
difficult  to  obtain  time -dependent  solutions  to  Eq.  (5),  either 
analytically  or  numerically.  We  therefore  seek  an  approximate  solution 
to  Eq.  (5),  based  on  the  observation  that  the  material  dispersion  in 
optical  waveguides  is  generally  very  small.  Our  previous  discussion 
regarding  the  balance  between  dispersion  and  nonlinearity  suggests  that 
the  nonlinear  effects  will  also  be  small  in  the  guide  case.  We  thus 
expect  a  negligible  effect  from  the  nonlinearity  on  the  mode  structure  of 
the  guide.  Physically,  the  spatial  trapping  of  the  beam  and  its 
localization  within  the  guide  is  dominated  by  the  transverse 
inhomogeneity.  If  nonlinearity  were  absent,  the  pulse  in  the  guide  would 
experience  broadening  due  to  dispersion.  The  correspondingly  small 
effect  induced  by  the  nonlinearity  is  capable  of  balancing  this 
longitudinal  dispersion.  Moreover,  it  seems  reasonable  on  general 
grounds  to  assume  that  transverse  inhomogeneity  has  only  a  weak  effect  on 
the  longitudinal  propagation  of  a  pulse  along  a  narrow  guide.  Rather, 


the  main  effect  of  the  transverse  inhomogeneity  is  to  determine  the  mode 
structure  of  the  linear  guide,  with  the  longitudinal  propagation 
characteristics  of  the  modes  subsequently  being  modified  by  the 
nonlinearity.  This  approach,  which  was  first  introduced  by  JT,  leads  to 
an  approximate  solution 


ACS,t)  =  <p(f)d(2,t) 


where  the  longitudinal  characteristics  of  the  pulse  are  found  by  taking 
the  transverse  average  of  Eq.  (5).  Taking  any  of  the  linear  mode 
solutions  for  A(?,t),  multiplying  Eq.  (5)  by  (f),  using  Eq.  (6)  for 
A(r,t)  and  integrating  over  the  cross  section  of  the  guide,  results  in 
the  following  wave  equation  for0(z,t); 


=  0 


(7) 


In  eq.  (7)  we  define  the  fuction  f  (z)  as 

and 

d,  =  <#)  if)  Vr  <j>  (f) 


The  electric  field  can  now  be  written  as 


& 


t) 


(8) 

(9) 

(10) 


(11) 


The  effect  of  the  waveguide  is  now  included  solely  by  means  of  the 
parameters  d^^,  6^  and  f(z). 
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We  now  consider  the  longitudinally  homogeneous  guide.  Here  f(z)  =f 
and  Eq.  (7)  supports  a  particular  solitary  solution  for  ^(z,t)  S 
0  (t-z/v),  where  ^  ^ 

(12) 


and 


with  the  conditions: 

The  significant  difference  of  our  result  (give  in  Eqs.  (14)  from  the 

homogeneous  medium  case,  is  the  factor  d^  above,  which  is  a  direct 
consequence  of  the  transverse  inhomogeneity,  i.e.,  the  waveguiding 
mechanism,  for  example,  we  consider  the  condition  for  the  "bright" 
soliton  solution,  given  in  Bq.  (13)  when  the  nonlinear  factor  is 
larger  than  zero.  The  requirement  ir>  0  implies  that 


(14a) 


(14b) 


(15) 


For  the  present  case,  in  which  d^  will  be  negative,  Eq.  (15)  may  be 
satisfied  even  when  k^"  ^0.  This  contrasts  with  previous  predictions 
and  with  the  one  dimensional  soliton  solutions  which  predict  that 
"bright"  solitons  will  exist  only  if  k^"^  0. 

We  point  out  that  Eq.  (7)  admits  a  "dark"  soliton  solution,  first  ^ 
realized  by  HT.  Here  a  solution  of  the  form  sech^[J('(t-z/v)  ^ 

exists  for  0^  ^  1.  For  the  particular  case  of  dark  solitons,  where 
0^  is  take  to  be  unity,  one  obtains 


0-  00  It  [V  -  '2lv)\ 


(16) 
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and  the  corresponding  conditions  given  in  ref.  10. 

For  plane  waves  k^"  =  0  indicates  the  zero  dispersion  condition  which 
determines  the  frequency  at  which  a  pulse  of  "^ero"  amplitude  will 
propagate  without  distortion.  In  our  case  the  zero-dispersion  condition 


becomes 


for  either  type  of  soliton.  Since  0  require  k^"  *>  0  for  the 

zero  -dispersion  condition  to  be  realized.  To  illustrate  these  results 
we  consider  the  truncated  quadratic  profile 

Cfp)  --  1  y  f  (II 

We  approximate  ^ (f^)  by  the  transverse  wave  solution  for  the  linear 
problem  when  ,  in  order  to  avoid  numerical  computations.  This 

approximation  is  very  good  for  the  lowest  modes.  Following  our 
procedures  described  above,  we  obtain  for  the  bright  soliton, 

0-  Oa  a  ~ 

The  zero-dispersion  condition  become$ 


-.-I 


The  zero-dispersion  condition  become  $ 


For  1  yum  light  in  fused  silica  we  then  have  k^L '^2  x  10^. 

If  f o  ~  example,  4.5yum,  while  if 

J’q  '^/\}  =  fo  3.2  yum.  These  estimates  become  increasingly 

inaccurate  as  r  /L  decreases.  JT  have  obtained  similar  results  in  a 
o 
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previous  paper,  where  regimes  of  bright  and  dark  solitons  as  a  function 

of  L  and  are  indicated  graphically  as  well.  JT  also  present 

calculations  for  step-index  fibers,  utilizing  the  transverse  function  of 

the  linear  problem  (i.  e.  the  Bessel  function).  We  may  use  their 

graphical  results  to  estimate  the  soliton  pulse  power  density  for  typical 

cases.  We  find,  e.g.,  for  the  lowest-order  mode  in  a  silica  fiber  with 

core  radius  5 ^m  and  carrier  wavelength  1 ^m,  that  the  power  density 
5  2 

P  /vlO  W/cm  for  a  pulsewidth  of  “t-  1  ns.  For  longitudinally 

inhomogeneous  waveguides  the  nonlinear  pulse  propagation  becomes  further 
complicated  when  the  guide  refractive  index  varies  as  a  function  of  z. 
Logitudinal  inhomogeneity  is,  in  fact,  quite  common  and  stems  from 
compositional  and  diametral  variations,  as  well  as  microbending.  Such 
variations  can  be  modelled,  at  least  approximately,  by  retaining  a 
longitudinal  dependence  in  the  effective  guide  refractive  index,  f=f (z) , 
in  Eq.  (7) .  Unfortunately,  there  are  no  known  methods  of  solving  Eq. 

(7)  in  this  instance.  On  the  other  hand,  methods  for  treating  slowly 
varying  variations  have  been  developed,  within  the  adiabatic 
approximation,  for  one-dimensional  nonlinear  propagation  in  the  absence 
of  frequency  dispersion  (CHEN,  H.H. ,  et  al.,  1976  and  1978).  We  are  thus 
led  to  consider  an  approximation  in  which  dispersion  and  longitudinal 


inhomogeneity  are  assumed  decoupled;  specifically,  we  assume 

_  21. 

i.e.,  the  dispersive  terms  are  replaced  by  their  values  in  the 


(20) 


homogeneous  limit.  Thus,  we  retain  interactions  between  transverse 
inhomogeneity,  dispersion  and  nonlinearity,  on  the  one  hand,  and  between 


I 


longitudinal  inhomogeneity  and  nonlinearity  on  the  other. 

The  solution  of  Eq.  (7)  in  this  instance  is  straighforward  but 
tedious.  The  reader  intersted  in  details  is  directed  to  BENDOW  et  al. 
(1979  and  1980) ;  we  will  here  just  state  the  results.  As  in  the 
longitudinally  homogeneous  case,  E  is  given  formally  by  Eq.  (11) ,  but 
now  takes  the  form  . 

Where  “Yk^k^'f/q,  and  are  time  dependent 

effective  wave  vector,  frequency  and  soliton  group  velocity,  of  the 
longitudinaly  inhomogeneous  waveguide,  respectively.  These  are  given  by 
(GIANINO,  P.D.,  et  al.,  1980) 

The  principal  result  of  the  above  development  is  that,  remarkably, 
the  soliton  retains  its  shape  in  the  presence  of  weak  longitudinal 
inhomogeneity,  although  its  velocity  change  with  time.  The 
time -dependent  velocity  in  this  case  contrasts  markedly  with  that  of  the 
longitudinally  homogeneous  case  where  v=q/k^kjj'f  is  a  constant.  For 
example,  when  the  inhomogeneity  is  nearly  linear,  i.e.,  when  An  =  2/2, 
then  !  ^ 

k«Ko  f 

so  that  the  soliton  acquires  a  constant  acceleration  proportional  to  •(.. 

2  2 

When  the  inhomogeneity  varies  quadratically , ^ n=  /  z  ,  then 


-X- 

k'oko 


"  IJt  / 


.1 


i.e.,  the  velocity  varies  sinusiodally  as  a  function  of  t.  Not 
surfSr isingly ,  the  soliton  executes  oscillatory  motion  as  a  function  of 
time,  characteristic  of  a  particle  trapped  in  a  harmonic  well.  Formal 
solutions  may  be  written  down  for  a  variety  of  other  cases  as  well,  such 
as  A  sin  JL  z,  in  which  case  v^  involves  elliptic  integrals 
(GIANINO,  P.D.,  et  al.,  1980).  In  this  instanced,  one  again  finds  that 
the  soliton  executes  a  complicated  oscillatory  motion  as  a  function  of 
time. 

REFERENCES 

1.  BENDOW,  B.  and  GIANINO,  P.D.,  1979  "Theory  of  Nonlinear  Pulse 
Propagation  in  Inhomogeneous  Waveguides",  Opt.  Lett.  j4,  164-166. 

2.  BENDOW, B.,  GIANINO,  P.D.,  TZOAR,  N.  and  JAIN,  M.  1980,  "Theory  of 
Nonlinear  Pulse  Propagation  in  Optical  Waveguides",  J.  Opt.  Soc.  Am.  70, 
539-546. 

3.  CHEN,  H.H.  and  LIU,  C.S.,  1976  "Solitons  in  Nonuniform  Media",  Phys. 
Rev.  Lett.  22,  693-697. 

4.  CHEN,  H.H.  and  LIU,  C.S.,  1978,  "Nonlinear  Wave  and  Soliton 
Propagation  in  Media  with  Arbitrary  Inhomogeneities",  Phys.  Fluids  21, 
377-380. 

5.  CHIAO,  R.Y. ,  GARMIRE,  E. ,  and  TOWNES,  C.H.,  1964,  "Self -Trapping  of 
Optical  Beams",  Phys.  Rev.  Lett.  12,  479-482. 

6.  CROSIGNANI,  B. ,  PAPAS,  C.H.  and  DI  PORTO,  P. ,  1980,  "Role  of 
Intensity  Fluctuations  in  Nonlinear  Pulse  Propagation",  Opt.  Lett.  5, 
467-468. 

7.  GERSTEN,  j.l.  and  TZOAR,  N.,  1975,  "Propagation  of  Localized 

Electromagnetic  Pulses  in  Plasma"  Phys.  Rev.  Lett.  35,  934-937. 


-10- 


8.  GIANINO,  P.D.  and  BENDOW,  B. ,  1980  "Simplified  Formulae  for  Solitons 

in  Media  with  slowly  Varying  Inhomogeneity",  Phys.  Fluids  220-221. 

9.  GLOGE,  D.,  1971,  "Dispersion  in  Weakly  Guided  Fibers",  ^pl.  Opt.  10 
2442-2445. 

10.  HASEGAWA,  A.  and  TAPPERT,  F. ,  1973,  "Transmission  of  Stationary 
Nonlinear  Optical  Pulses  in  Dispersive  Dielectric  Fibers  I.  Anomalous 
Dispersion",  ^pl.  Phys.  Lett.  22,  142-144.  "II.  Normal  Dispersion", 
appl.  Phys.  Lett.  22'  171-172. 

11.  JAIN,  M.  and  TZOAR,  N. ,  1978,  "Propagation  of  Nonlinear  Optical 
Pulses  in  Inhomogeneous  Media",  J.  Appl.  Phys.  4_9i  4649-4654,  and, 
"Nonlinear  Pulse  Propagation  in  Optical  Fibers",  Opt.  Lett.  2»  202-204. 

12.  KAW,  P.K.,  NISHIKAWA,  K.,  YOSHIDA,  Y.  and  HASEGAWA,  A.,  1975, 
"Two-Dimensional  and  Three-Dimensional  Envelope  Solitons",  Phys.  Rev. 
Lett.  25,  88-91. 

13.  KECK,  D.,  1976,  in  "Fundamentals  of  Optical  Fibers  Communications", 
edited  by  M.  Barnoski,  Academic  Press. 

14.  SCOTT,  A.C.,  CHU,  F.  Y.  F.  and  MC  LAUGHLIN,  D.W. ,  1973,  "The  Soliton 
A  New  Concept  in  Applied  Science",  Proc.  IEEE  61,  1443-1483.  An 
extensive  list  of  references  in  contained  in  this  review  article. 

15.  SHIMIZU,  F. ,  1967,  "Frequency  Broadening  in  Liquids  by  a  Short  Light 
Pulse",  Phys.  Rev.  Lett.  19,  1997-1100. 

16.  TZOAR,  N.  and  GERSTEN,  J.I.,  1976,  in  "Optical  Properties  of  Highly 
Transparent  Solids",  edited  by  S.S.  Mitra  and  B.  Bendow,  Plenum  Press. 

17.  WILCOX,  J.Z.  and  WILCOX,  T.J. ,  1975,  "Stability  of  Localized  Plasma 
Model  in  Two  and  Three  Dimensions",  Phys.  Rev.  Lett.  34,  1160-1163. 

18.  ZAKHROV,  V.E.  and  KUZNETSOV.  E.  A.,  1979,  "Three-Dimensional 

Solitons",  Sov.  Phys.  JETP  21'  285-286. 


-11- 


SELF-PHASE  MODULATION  IN  OPTICAL  WAVEGUIDES 


Self-phase  modulation  was  first  observed  by  Shimizu^  when  a 

modulated  spectrum  appeared  after  self -focusing  had  taken  place  in  a 

liquid-filled  cell  and  was  explained  as  phase  modulation  due  to  the 

2 

intensity-dependent  refractive  index.  It  has  since  been  observed  in 
the  absence  of  self -focusing  or  self -trapping  and  at  low  powers  by  using 
liquid-filled  glass  fibers.  Recently,  some  measurements  of  frequency 
broadening  of  mode-locked  laser  pulses  due  to  self-phase  modulation  in 
single-mode  silicacore  fibers  have  been  reported.^  The  theory  used  in 
ref.  3  is  the  one  developed  first  by  Shimizu  for  short  samples,  and 
shows  symmetic  broadening  of  the  spectia  due  to  self-phase  modulation. 

We  present  calculations  for  the  frequency  broadening  due  to 
self-phase  modulation  which  show  an  asymmetry  in  the  output  spectrum  for 
symmetric  initial  pulse.  Comparison  with  results  of  Ref.  3  is  difficult 
to  make  at  present,  firstly  due  to  large  asymmetry  of  the  incoming  pulse 
and  secondly  due  to  dispersion  effects  which  will  be  discussed  later  in 
this  paper.  We  thus  limit  ourselves  to  the  development  of  the  theory  and 
determine  the  realistic  conditions  undr  which  comparison  between  theory 
and  observation  is  easy  to  make. 

Wo  limit  ourselves  to  a  one-dimensional  wave  propagation  of  an 
optical  pulse  propagating  in  a  glass  characterised  by  a  nonlinear 


refractive  index  n.  The  electric  field  E  (  z,t)  is  given  by  the  wave 
4  5 

equation  ' 

}‘E  _i  2-E  =  fiBi'-e)  , 


-12- 


*1 


where  n2^  represednts  the  nonlinear  part  of  the  refractive  index 

K(uj,e)  -  ' 

with  =n  ( )  ,  where  (a)^  is  the  frequency  of  the  electric  field  and  is 
the  linear  displacement. 

We  write  the  electric  field  in  terms  of  a  slowly  i^arying  envelope 

e  i  in 


where  q  is  the  propagation  constant  and  find  that  A(2,t)  obeys  the 
equation  h 


where 


In  Eq.  (3)  we  retain  only  the  first  and  second  time  derivatives  of 
A.  Here  we  make  the  usual  assumption  that  our  dielectric  is  weakly 
dispersive.  This  is  an  extremely  good  approximation  for  Si02  glasses 
for  our  operating  frequency  which  is  much  smaller  than  the  electronic 
resonance  frequency  and  much  larger  than  the  ionic  resonance  frequency  of 


S1O2  . 

To  first  order  we  take  q=kQ  as  for  the  plane-wave  situation,  and 
identify  k^'  to  be  the  reciprocal  of  the  group  velocity  Vg,Thus  we 
are  left  with  the  following  equation; 
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In  treating  self-phase  modulation  one  usually  replaces  the  right-hand 
side  of  Eq,  (5)  by  -2(n2  /no)k^  lA/'A  as  the  dominant  contribution, 
assuming  that  Here  T  represents  the  width  of  the  pulse. 

However,  as  we  shall  see  shortly,  this  assumption  is  correct  only  for 
short  enough  samples.  For  long  samples,  as  for  the  case  of  self-phase 
modulation  in  optical  waveguides,  we  have  to  consider  the  corrections 
arising  from  the  time  derivatives  of  /  a|  A  .  We  will  treat  this  effect 


perturbatively  and  thus  will  retain  only  the  first-order  term,  i.e.,  the 
first  derivative  of  /  a|  A  .  Equation  (5)  is  therefore  approximated  by 


We  now  make  a  coordinate  transformation  to  a  moving  coordinate  system 
defined  by 

^  =  2  ; 


Equation  (6)  transforms  to  yield 


To  further  simplify  Eq.  (7)^ we  realized  that  without  the  dispersion 

n 

(k^)  and  the  nonlinearity  (02)  the  solution  for  A  is  given  by  an 

arbitrary  function  ofZ^,  say,  F[  (z  -  v  t)  /z“  ] ,  where  ?  is  the  legth 

^  ! 

of  the  pulse.  It  then  becomes  obvious  that  while  and  kg  P/p-j.  are 

of  the  order  of  ‘z“^  ,  the  term  kg  is  proportional  to  We  can  now 
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use  the  slowly  changing  envelope  approximation  which  required  "z  to  be 
much  larger  than  so  that  (  J -  2  ^ )  A  is  neglected 

relative  to  ("o  ^i.e.,  the  effects  of  back-scattered  radiation  are 

neglected.  This  implies  that  the  changes  in  A  per  wavelength  are 
extremely  small.  This  condition  is  compatible  with  our  suggested 
experimental  situation,  where  changes  in  A  are  only  observed  after  the 
pulse  propagates  hundreds  of  meters  in  the  guide.  We  next  define  the 
constants 


If 


and  obtain  our  final  equation 


(81 


(91 


Dt' 


The  detailed  solution  for  A  is  given  in  the  paper  by  Tzoar  et  all  (Phys 
Rev  A  23  1266,  1981) .  We  present  here  only  the  approximate  result  for 
the  frequency  spectrum  of  the  electric  field  having  a  Gaussian  envelope. 


E  ( =  e  A .  e‘  j"p 


(10) 


where 


F=expi 


(Ill 


In  the  above  expansion,  terms  proportional  to  d.  arise  from  dispersion 
effects  and  are  seen  to  give  rise  to  a  symmetric  modulation  of  the 
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phase.  The  term  proportional  to  0  ,  however,  is  proportional  to"C-  and 
thus  results  in  asymmetric  phase  modulation. 

To  analyze  the  expression  in  Eq.  (18)  we  note  that  in  our  problem  we 
have  two  fundamental  lengths,  the  nonlinear  length  (  ^Ao)  and  the 
dispersive  length  5 j =T/  Ml  ,  and  also  the  dimensionaless  parameter  (Oq  T. 
We  rewrite  Eq.  (11)  in  terms  of  these  parameters  as 


Here  =t3^T^/16  and  L2  effect  of  the  linear 

dispersion  can  be  omitted  in  many  cases  since  it  is  easy  to  rnake^^^?^ 
The  correction  terms  depend  on  and  L2  .  The  term  is  the 
correction  arising  from  the  finite  duration  of  the  pulse  and  is  the  only 
term  that  contributes  in  our  theory  to  asymmetric  spectrum.  To  observe 
the  asymmetric  effects,  we  want  the  term  to  be  the  dominant 
correction.  We  thus  need  to  be  larger  than  . 

We  note  that  the  ratio  /L2  should  be,  ideally,  smaller  than 

unity  for  our  purposes.  Here 


is  independent  of  T  and  becomes  small  either  at  high  beam  intensity  or 

•f 

for  vanishingly  small  dispersion,  i.e.,  k^CrO. 

The  dispersion  effects  can  be  reduced  by  operating  at  a  frequency 
where  the  dispersion  is  negligible^  as  has  been  recently  demonstrated. 
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Thus,  in  Eq.  (17) ,  the  integrand  reduces  to 


F  =  1 


^  AoC 


(13) 


In  order  to  make  an  estimate  of  the  asymmetric  effect,  we  calculate  the 

value  of  the  coefficient  8  ■fix  Aq/T^  ,  using  values  that  can  be 

Ns  13 

attained  in  the  laboratory.  Taking  5=1  km,  X'^T=5ps,  =1.4  10 

7  ^22 

esu,  and  Aq  =500  S  V/cm,  we  calculate  S^iXAq/T  /v  0.25. 

Thus  the  effect  of  the  asymmetric  term  compared  to  1  is  about  25%  for  the 
parameters  chosen  above  and  should  be  experimentally  measurable. 

We  would  like  to  point  out  that  when  information  transmision  depends  on 
the  phase  of  the  carrier,  self-phase  modulation  must  be  considered  as  a 
noise  source  that  may  bleach  the  signal.  In  this  work  we  have  developed 
a  theory  which  is  better  suited  for  long  (^1  km)  optical  waveguides.  We 
have  shown  that  self-phase  modulation  will  result  in  asymmetric  power, 
and  similaly,  in  asymmetric  intensity  spectra. 
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7.  For  a  waveguide  with  a  cross  section  of  10  cm  the  peak  power 
is  about  2.5  W. 
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CONTINUOUS-WAVT';  PROPAGATION  IN  NONLINEAR  OPTICAL  WAVEGUIDES 


1.  INTRODUCTION 

Transmission  of  intense  light  beams  through  waveguides  is  of  interest  in  a 
wide  variety  of  potential  applications,  such  as  integiated  optics  (four-wave 
mixing,  optical  bistable  devices),  medical  procedures  (surgery, 
cauterization),  industrial  processing  (cutting,  welding),  and  power 
transmission  (for  electrical  hazard  zones,  say) .  The  nonlinearity  of  the 
refractive  index  of  optical  waveguides  affects  both  the  spatial  and  the 
temporal  characteristics  of  intense  beams  propagating  along  the  guide.  Here 
we  investigate  certain  spatial  characteristics  of  cw  beams  in  waveguides 
possessing  a  real  nonlinear  dielectric  constant.  In  other  work  various 

temporal  effects  in  nonlinear  waveguides,  including  soliton  propagation^ 

2 

and  self-phase  modulation,  have  been  analyzed.  Various  phenomena 
associated  with  the  imaginary  part  of  the  nonlinear  dielectric  constant  have 
been  discussed  previously.^ 

The  inhomogenous ,  nonlinear  wave  equation  describing  beam  propagation  in 
waveguides  is  not,  in  general,  amenable  to  analytic  solutions,  nor  is  it 
especially  well  suited  to  numerical  methods  of  solution.  For  these  reasons, 
we  have  chosen  to  introduce  several  simplifying  approximations  commonly 
employed  in  the  literature  that  make  the  wave  equation  analytically 
tractable.  Specifically,  we  restrict  consideration  to  Gaussian  beams  and 

A 

employ  the  widely  used  parabolic  or  paraxial  approximation.  This  is 
essentially  a  ray  treatment,  which  is  most  accurate  in  the  vicinity  of  the 
beam  axis.  Such  an  approach  has  been  employed  extensively  in  previous  work 
to  investigate  nonlinear  focusing  in  homogenous  media, as  well  as  linear 
propagtion  in  inhomogeneous  media. 
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In  the  present  investigation  we  extend  the  paraxial  approximation  to 
encompass  simultaneously  inhomogeneity  as  well  as  nonlinearity,  but 
restrict  ourselves  to  the  effects  of  nonlinearity  on  the  focal 
characteristics  and  on  the  mode  mixing  of  Gaussian  beams  that  are 
incident  on-axis.  Use  of  these  approximations  and  simplifications  allows 
us  to  examine  some  of  the  principal  physical  effects  associated  with 
nonlinear  beam  propagation  and  helps  to  reduce  mathematical  complexity  to 
a  minimum. 

2.  FOCAL  CHARACTERISTICS  OF  GAUSSIAN  BEAMS 

Here  we  derive  an  expression  for  the  electric  field  of  a  cw  Gaussian  beam 
in  a  nonlinear  waveguide  by  utilizing  the  paraxial  approximation.  The 
results  obtained  here  also  incorporate  guide  attenuation  and  longitudinal 
dielectric  constant  variations.  Since  it  is  not  possible  to  solve  the 
governing  equations  analytically  in  the  general  case,  we  subsequently 
choose  to  specialize  in  the  case  of  lossless,  longitudinally  uniform 
guides.  In  this  instance,  one  is  able  to  obtain  simple  analytic 
solutions  that  clearly  reveal  the  effects  of  waveguiding  and 
nonlinearity  on  cw  beam  propagation.  We  developed  our  theory  for  an 
azimuthally  symmetric  waveguide  with  the  z  axis  along  the  length  of  the 
guide  and  the  radial  coordinate  in  the  transverse  direction.  We  choose 
a  model  dielectric  constant  of  the  form 

6(f,z,£)-=  £„W- 

In  Eq.  (1),  waveguiding  is  introduced  through  the  quadratic  term  in^^. 
This  widely  used  approximation  assumes  that  the  beams  under  consideration 


are  tightly  bound  within  the  core  of  the  guide.  The  nonlinearity  is 


introduced  through  a  term  in  6  proportional  to  the  square  of  the 
electric  field  E,  i.e.,  to  the  field  intensity. The  imaginary  term 
accounts  for  guide  attenuation  (from  absorption  and  scattering,  say). 
Finally,  we  have  allowed  for  longitudinal  variations  in  the  dielectric 
constant  by  incorporating  a  z  dependence  into  and 

Thus,  despite  the  simplicity  of  this  model  for  the  dielectric 
constant,  one  can  account,  at  least  approximately,  for  any  combination  of 
inhomogeneity  (both  transverse  and  longitudinal) ,  nonlinearity,  and 
attenuation. 

If  we  assume  a  time-harmonic  dependence  of  the  electric  field,  E 

;u,t 

e  ,  then  the  time-independent  scalar  wave  equation,  neglecting  gradient 
terms,  takes  the  form^  (k^  =  OJ/c) 

[  ^  .  (2) 
We  look  for  azimuthally  symmetric  solutions  of  the  form 


Where  k^(z)  =  ^q{z)  .  By  substituting  Eq.  (3)  into  Eq.  (2) 

and  neglecting^  the  terms  proportional  to  C'O 

OmJ  JL  •  O"®  obtains 

(d.iK+2.ik)^  =  i- 

Equation  (4)  constitutes  the  starting  point  for  the  paraxial 
approximation,  in  which  variations  in  A  as  a  function  of  z  are  assumed  to 
be  negligible  compared  with  a  wavelength.  It  reduces  to  the  same 
starting  equation  employed  by  various  other  authors^*^  when  , 

ki, 

and/or  £  are  set  equal  to  zero. 
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In  ocder  to  proceed  with  the  solution  of  the  wave  equation,  one 
further  expresses  A  in  terms  of  a  real  amplitude  Ag  and  phase  S  as 


(5) 


and  substitutes  Eq.  (5)  into  Eq.  (4)  One  then  obtains  two  simultaneous 
equations  for  A.  and  S; 


2  S 
k 


Si + ^  IT  +(-jf j  =  T?A.(  rp "  f 

These  equations  become  soluble  if  we  employ  for  the  nonlinear  term  the 
quadratic  approximation  introduced  by  Ref.  8;  namely,  we  take 

iAi^  r’' 

lAl  =  lAl^  t  J 

Where  the  subscript  0  indicates  evaluation  at  ^  =0.  One  then  finds,  in 
analogy  with  Refs.  4-7,  that  Eq.  (6)  admits  the  Gaussian  beam  solution 
(kj  =ik  6j/£p),  ^ 

(8b) 

Where  f  and  ^  satisfy 

4X  +  -J—  ^  =  (r;'- 

^  ^  h'L  (9a) 


and 


(10) 


A»<,.Eo'1[  WA.(o)1 

4 

If  the  incident  beam  is  a  plane  wave  /  then  the  boundary  conditions 
accomplanying  Eqs.  (9a)  and  (9b)  are 

The  complete  solution  for  E  is  then 

E  =  -It, 

n  (12) 

Where  f  and  (|)  are  determined  from  Eqs.  (9)  and  (11).  This  set  of 
equations  encompasses  a  complete  formal  solution,  within  the  paraxial 
approximation,  for  the  electric  field  of  cw  Gaussian  beams  in  a  nonlinear 
wavequide  described  by  the  model  dielectric  constant  of  Bq.  (1) . 

Certain  general  characteristics  of  the  propagation  are  evident 
directly  from  inspection  of  Eqs  (9) -(12).  For  example,  the  spatial 
profile  of  the  beam  is  observed  to  be  a  Gaussian  profile  whose  focal 
parameter  varies  with  distance  along  the  guide  z.  The  z  dependence  of 
the  focal  parameter  is  influenced  by  the  transverse  and  longitudinal 
variations  of  the  dielectric  constant  of  the  guide,  the  nonlinearity,  and 
the  loss.  Equation  (12)  indicates  that,  oecause  of  the  loss  term  in  G , 
the  field  amplitude  decreases  exponentially  with  propagation  distance  in 
a  familiar  Beers  law  manner.  Finally,  note  that  the  phase  of  the  wave 
acquires  a  somewhat  complicated  z  dependence. 

The  specific  characteristics  of  the  propagation  are  obtained  once  the 
differential  equations  for  f  and  S  have  been  solved.  Usually,  the 
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primary  quantity  of  interest  is  the  beam  intensity,  in  which  case  phase 
factors  are  irrelevant  and  only  the  solution  for  f(z)  is  required. 

Clearly,  the  solution  of  the  nonlinear  inhomogeneous  differential 
equation  for  f  [Eq.  (9a)]  cannot,  in  general,  be  obtained  by  analytic 
means  but  requires  the  use  of  tedious  numerical  techniques  instead.  For 
this  reason,  we  choose  not  to  treat  the  general  case  here;  readers 
interested  in  the  application  of  numerical  techniques  to  equations  of  the 
form  of  Eq.  (9a)  are  directed  to  Ref.  7.  Rather,  in  the  remainder  of 
this  report,  we  specialize  in  the  simplified  limit  of  a  lossless, 
longitudinally  uniform  waveguide.  In  this  limit,  is  is  possible  to 
obtain  explicit  analytic  solutions  for  f(z)  that  clearly  reveal  the 
effects  of  waveguiding  and  nonlinearity  on  cw  beam  propagation. 

Although  in  the  following  we  especialize  to  lossless,  longitudinally 
uniform  guides,  it  should  be  pointed  out  that  Eqs.  (9) -(12)  provide  a 
convenient  starting  point  for  treating  more  general  cases.  The  solution 
of  the  wave  equation  (2)  in  the  simultaneous  presence  of  nonlinearity, 
transverse  and  longitudinal  inhomogeneity,  and  loss  is  prohibitively 
difficult.  If  the  paraxial  approximation  is  employed,  the  problem  is 
reduced,  essentially,  to  the  solution  of  a  single  second-order  ordinary 
differential  equation  [Eq.  (9a)). 

Solutions  for  Longitudinally  Uniform,  Lossless  Waveguides: 

For  a  longitudinally  uniform  lossless  guide,  we  set^(z)  r  ^  r,  (z)= 

O  OR  ^ 

,  and  6^=0.  Then  the  focal  parameter  f  (z)  satisfies  a  simplified  equation 
of  the  form 
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Here  represents  an  effective  radius  determined  by  the  relaive  sizes 
of  diffraction  and  nonlinearity. 

2 

Consider  first  the  case  in  which  0  (R^  ^ 


solution  for  f  yields: 


C=6o/£R?e  y  T-  'Z. 


(14) 


Thus,  for  sufficiently  small  nonlinearity,  the  beam  remains  "trapped"  in 
the  waveguide  and  continues  to  oscillate  with  spatial  per  iod^^  ( , 
despite  the  nonlinearity  of  the  medium.  At  this  level  of  approximation, 
then,  the  guide  effect  totally  dominates  the  diffraction  effect,  whereas 
nonlinearity  influences  only  the  amplitude  of  the  oscillation  but  not  its 
period . 


Several  observations  are  useful:  First,  it  is  easily  shown  that 
for  the  above  cases  (R^'7  0) ,  f  never  vanishes,  i.e.,  the  beam  never 
collapses.  Second,  we  note  from  Eq.  (21)  that,  if  C=l,  then  f=l 
independently  of  z.  From  Eqs.  (13)  and  (21)  this  occurs  when 


Here  Wj^  is,  in  fact,  just  the  spot  size  of  the  lowest-order  eigenmode 
of  an  approximate  quasilinear  wave  equation  [one  obtained  by 
approximating,  the  nonlinear  term  by  using  Eq.  (7)  at  z=0].  Thus, 
within  this  approximation,  if  the  spot  size  of  the  incident  beam  matches 
that  of  tfhe  lowest-order  mode  of  the  guide  at  z=0,  then  one  predicts 
(not  surprisingly)  that  the  beam  will  propagate  unchanged  as  a  function 
of  z. 

2 

Now  consider  the  case  in  which  Rg  ^0  (equivalently,  C<0).  In 
this  instance  we  obtain  a  singular,  self -focusing -type  solution  for  f, 
similar  to  the  corresponding  homogeneous  nonlinear  medium  case 
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irrespective 
find  that 


f 


of 

Z 


the  waveguide  effect.  In  particular,  in  this  case  we 

ih-  ici)-+  0+ 


Here  we  define  a  self -focusing  length  (Zpj^in  dimensionless 

2. 

units)  determined  by  the  condition  f  C^)  =0,  whence 


(17) 


2 

One  obtains  a  periodic  solution  for  f  ,  which  may  be  extended  to  the 
region  3  V  as 

^  (^/-f /C/)  (18) 

At  and  "S  =  3  "Tp  we  obtain  f  =0,  and  the  derivative  of  f  is 

2  2 
discontinuous.  The  solution  for  f  may  be  summarized  as  follows:  f 

is  given  by  Eq.  (16  )  for  tf  ,  and 

2 

for  j=l  1,3,5...  This  solution  indicates  that  f  =  1  for  even  multiples 
_  2 

of  3p  and  that  f  =0  (beam  collapse)  for  odd  nultiples  ofjp*  In  this 
case,  nonlinearity  dominates  both  waveguiding  and  diffraction,  resulting 
in  self -focusing . 

Just  as  in  the  homogeneous  medium  case,  the  focusing  of  the  beam  to  a 
point  (f-^  0)  is,  in  the  waveguide  case,  a  consequence  of  the  insta¬ 
bility  induced  by  nonlinearity.  The  critical  power  for  self -focusing 
[see  Eq.  (27)  below],  as  determined  from  the  condition  R 
oscillatory  nature  of  the  solution  both  remain  unaltered  by  the 
waveguide.  The  main  effect  of  the  guide  is  to  modify  the  functional  form 
of  f  and  alter  the  focal  length  Zp  .  But,  overall,  effects  that  are 
due  to  nonlinearity  dominate  over  those  that  are  due  to  waveguiding. 


and  the 


26- 


In  order  to  estimate  the  power  level  P  at  which  the  influence  of 

2 

nonlinearity  becomes  significant,  we  first  reexpress  E*  in  esu  as 

£■/  =  8  ?/^ 

As  in  the  nonguided  case,^  we  can  consider  the  critical  power  P  for 
which  nonlinearity  exactly  balances  diffraction.  This  occurs  when  C=0, 
i.e.,  the  f  ^  term  in  Eq.  (13)  vanishes  identically,  whence  the  beam 
propagation  is  determined  exclusively  by  waveguiding  (depending  only 
on  /(C  )  .  Setting  C=  0  yields  (in  esu)  for  P^  : 

which  is  independent  of  the  initial  beam  radius.  By  using  Eqs.  (27)  and 

(28),  it  is  possible  to  reexpress  the  condition  for  waveguiding  as 

P^P  whereas  the  condition  for  self-focusing  becomes  P^P  . 
c  c 

Clearly,  the  influence  of  nonlinearity  on  the  focal  parameter  is  greatest 
for  with  only  small  changes  occurring  for  P<Jj'P^. 

The  actual  value  of  the  critical  power  P^  varies  considerably  from 

case  to  case,  depending  on  the  propagation  medium  and  the  operating 

\  11 
wavelength.  For  for  example,  P^  '^lO  esu 

[10  kW]  for  materials  with  a  large  third-order  nonlinear  susceptibility 

(various  semiconductors  and  liquids,  such  as  CS2) ;  on  the  other  hand, 

P^  ^  10^^  esu  [1  MW]  for  weakly  nonlin‘'^.r  materials  (such  as 

silicate  glasses) . 
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3.  NONLINEARITY-INDUCED  MODE  MIXING 

Consider  the  case  of  a  cw  Gaussian  beam  propagating  in  a  lossless, 
longitudinally  uniform  waveguided.  We  found  that,  for  beam  powers  below 
the  critical  power,  the  oscillation  period  of  the  focal  parameter  remains 
unchanged,  but  its  amplitude  becomes  a  function  of  the  incident  beam 
power.  In  applications  in  which  the  spatial  mode  structure  of  beams  is 
significant,  such  as  phase  conjugation  of  images  in  fiber  guides^^  by 
four-wave  mixing,  it  is  useful  to  determine  the  mode  mixing  (i.e., 
variation  in  the  modal  content  of  the  beam  as  a  function  of  z)  that 
accompanies  the  dependence  of  f  on  beam  power.  A  description  in  terms  of 
mode  mixing  makes  the  most  sense  below  the  critical  power,  where  stable 
waveguided  propagation  occurs  (although  one  may  formally  calculate  the 
mode  mixing  accompanying  self-focusing  as  well). 

We  choose  for  consistency  to  define  mode  mixing  in  terms  of  an 
approximate  initial  set  of  modes  obtained  within  the  spirit  of  the 
paraxial  formulation. Specif ically,  if  we  utilize  Eq.  (7)  and  (8)  to 
evaluate  the  nonlinear  term,  then  we  obtain  what  will  be  referred  to  as  a 


quasi-nonlinear  wave  equation  for  the  transverse  eigenfunction 
and  corresponding  propagation  constants  k  at  the  fixed  point  z=0: 

(k*  -k'juWf) 

where  we  have  used  f(0)  =1.  The  eigensolutions  of  this  equation  are^ 


(wr iting  for  &p{0)  and  (r^^for  <£p^( 0)  ) 

w,'%  1/ 6y,  Eo  lu  f!'  ;  w;^  k  / 6. 
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(24) 


where  L  is  the  Laquetre  polynomial  and  n  takes  on  integral  values, 

V 


i<‘=  kV'tfiw  £«/<£.;- "I 


For  tne  limit  of  a  linear  guide  on  simply  sets^^O,  in  which  case  Eqs. 
(3)  and  (31)  reduce  to  the  well-known  expressions  for  azimuthally 

symmetric  transverse  modes  of  a  linear  longitudinally  homogeneous 

.  ,  10 
guide. 

In  analogy  to  Ref.  (4),  we  expand  the  expression  for  E/Eq  in  Eg. 
(12)  in  terms  of  the  u  's  as 


(25) 


where  ^^(j*)  is  given  by  Eg.  (23)  Employing  the  orthonormality  of  the 
“U^'s  yields  after  much  algebra  the  desirable  result. 

By  defining  a  dimensionless  mixing  coefficient  as: 


one  obtains  for  the  guasilinear  expansion  set 


Where : 


.  i  •  r  - 


(26) 


(27) 


(28) 


Several  general  observations  regarding  mode  mixing  in  nonlinear 
waveguides  follow  directly  from  an  inspection  of  Eg.  (37).  For  example, 
we  see  that  mode  mixing  is  generally  present.  Mode  mixing  is  absent  only 
in  the  special  case  in  which  C  »1  (equivalently,  0^  +  ()^  =1),  in  which 
instance  f (z)  =  constant.  Physically,  this  represents  the  case  in  which 
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the  incident  beam  spot  size  exactly  matches  that  of  the  lowest-order 
eigenmode  of  the  guide  at  z=0  . 

For  fixed  (Tj  and  ,  the  mode  mixing  is  observed  to  increase  with 
increasing  mode  number  n.  The  dependence  on  0^  and  0^  is  somewhat  more 

complicated.  Clearly,  for  large  values  of  7',  ) ,  mode  mixing 

increases  as  the  nonlinearity  increases  (i.e.,  as  P  P  ).  However, 
another  determining  factor,  as  discussed  above,  is  the  proximity  of  (/f 
to  unity.  When  (r^  +  0^=l,  mode  mixing  vanishes,  but,  as  the  parameters 
depart  from  this  condition,  the  mode  mixing  increase. 

As  is  the  case  for  P^  ,  the  parameters(J^and  fj^can  take  on  a  wider 
range  of  values  depending  on  the  operating  wavelegth,  the  propagation 
medium,  and  the  geometrical  parameters  of  the  beam  and  waveguide.  To 

obtain  some  indication  of  typical  values,  let  us  again  consider  ^^=^m 

and  utilize  the  range  of  values  indicated  previously  for  P 

o 

(10^-10^  W) .  Then 

<r,  ~  (/o'‘'-/o'OF 

where  P  is  the  number  of  watts. 

Typical  ranges  of  values  for  the  waveguide  and  beam  parameters 
are  ^  -10  ^  and  (  |<^ )  ^-10 respectively. 

The  corresponding  range  for0^is  then  10”^-10^^,  although  the 

lower  portion  of  this  range  is  more  typical  of  standard  fiber  guides. 

4.  SUMMARY  AND  DISCUSSION 

In  this  section  we  have  examined  various  aspects  of  propagation  of  cw  beam 
in  nonlinear  optical  waveguides  and  contrasted  the  results  with  those  for 
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linear  waveguides  and  homogeneous  nonlinear  media.  Analytic  expressions 
for  the  dependence  of  the  focal  parameter  f  of  a  Gaussian  beam  as  a 
function  of  distance  z  are  obtained  within  the  paraxial  approximation. 

We  have  demonstrated  the  existence  of  two  distinct  regimes.  For  powers 
below  the  critical  power  for  self -focusing  in  a  homogeneous  medium,  the 
waveguide  effect  dominates;  the  beam  becomes  trapped  and  merely  varies 
sinusoidally  with  z.  Above  the  critical  power,  nonlinearity  dominates 
the  propagation;  the  beam  becomes  unstable,  and  oscillatory  self -focusing 
behavior  is  predicted, in  strict  analogy  to  a  homogeneous  medium.  The 
detailed  solutions  reveal  that  in  the  waveguide-dominated  case,  the 
functional  form  for  f  and  the  self -focusing  distance  Zp  are  modified 
because  of  the  waveguide.  One  must  carry  the  approximations  to  higher 
order  in  order  to  obtain  changes  in  oscillation  period  or  critical  power 
because  of  interactions  between  waveguiding  and  nonlinearity. 

We  have  defined  mode  mixing  in  terms  of  the  quasilinear  modes  (which 
have  been  corrected  for  nonlinearity  in  an  approximate  fashion) .  We 
find  that  mode  mixing  is  always  present  in  nonlinear  guides  (except  for 
the  trivial  case  in  which  the  beam  focal  parameter  is  a  constant). 

These  results  contrast  strongly  with  those  for  the  linear  guide,  in  which 
mode  mixing  occurs  only  in  longitudinally  inhomogeneous  guides.  We 
obtain  analytic  expressions  for  the  mixing  coefficient  g(n,z)  for 
longitudinally  homogeneous  guides.  The  results  indicate,  in  general, 
increased  mixing  with  increasing  mode  number  and  with  increasing  beam 
powe  r . 
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Certain  aspects  of  the  present  work  may  be  useful  in  connection  with 
other  nonlinear  problems.  For  example,  consider  the  calculation  of 
nonlinear  pulse  propagation  in  waveguides.  All  treatments  of  the 
time -dependent  case  to  date  have  essentially  decoupled  the  longitudinal 
and  transverse  motions  by  averaging  over  the  transverse  coordinates  of 
the  guide. The  present  work  suggests  that,  for  appropriate  choices  of  the 
guide  parameters  (e.g.  for  sufficiently  low  beam  power),  the  transverse 
motion  is  indeed  strongly  waveguide  dominated,  with  only  negligible 
modifications  because  of  nonlinearity.  Moreover,  the  mixing  of  modes  of 
sufficiently  low  order  may  be  small,  in  which  case  the  concept  of 
individual  solitons  associated  with  the  separate  modes  of  a  multimode 
guide  may  be  valid  to  some  level  of  approximation.  In  any  case,  the 
present  formulation  may  help  assess  under  what  circumstances  it  is 
reasonable  to  decouple  transverse  and  longitudinal  effects  in  nonlinear 
pulse  propagation. 

Another  promblem  for  which  the  present  results  may  be  useful  is  that 
of  multiple-wave  mixing  in  waveguides.  Most  treatments  to  date  assume 
that  the  pump  beams  are  propagated  in  the  lowest-order  transverse  mode  of 
the  guide.  If  these  beams  do,  in  fact,  undergo  mode  mixing,  then  the 
spatial  characteristics  of  the  signal  beams  may  become  altered  as  well. 
The  present  work  provides  useful  input  regarding  the  extent  of  mode 
mixing  as  a  function  of  the  relevant  system  parameters. 

We  comment  briefly  on  the  possibility  of  experimentally  observing  the 
phenomena  discussed  in  this  paper.  In  this  connection,  one  must  dis¬ 
tinguish  between  the  waveguiding  (P<P_)  and  self-focusing  (  P>P  ) 
regimes.  The  relatively  small  changes  induced  by  nonlinearity  for  P^ 
are  significant  primarily  with  respect  to  mode  mixing.  They  will 
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therefore  be  evident  only  in  processes  that  are  highly  sensitive  to  mode 
structure,  such  as  modulation  of  spatial  signals  by  multiple-wave  mixing, 
proximity  coupling,  or  guided-wave  microbending  loss  sensing.  On  the 
other  hand,  the  self-focusing  that  is  predicted  for  will  be  a 

large  effect,  provided  that  it  is  not  overshadowed  by  other  competing 
nonlinear  processes.  Certain  related  nonlinear  phenomena  can  easily  be 
excluded  or  distinguished  by  their  special  characteristics.  For 
example, soliton  propagation  requires  pulsed  as  opposed  to  cw  operation 
and,  in  any  case,  is  a  relatively  low-power  phenomenon  that  serves  to 
balance  the  usually  smalll  dispersion  encountered  in  typical  waveguides. 
On  the  other  hand,  various  parametric  processes,  such  as  stimulated  Raman 
and  Btillouin  emission  and  multiple-wave  mixing,  can,  in  fact,  compete  or 
coexist  with  self -focusing,  depending  on  the  operating  frequency, 
propagation  medium,  and  other  system  parameters.  In  most  instances, 
such  processes  can  be  distinguished  for  self -focusing  because  they 
possess  different  dependences  on  propagation  distance.  The  growth  of 
parametric  processes  is  usually  assumed  to  vary  as  exp(gPz/A),  where  g  is 
a  gain  coefficient  and  A  is  the  cross-sectional  area  of  the  light  beam. 
Neglecting  all  other  factors,  let  us  define  a  minimum  distance 
for  parametric  processes  be  setting  gPZp  /A^l.  For  self -focusing ,  on 
the  other  hand,  when  P»  P^,  and  ^  6^  one  can  show  that 
/P)  .  Thus  we  find  that 

r  C 


Consider,  for  example,  fused  silica  at^*»l yim  and  P/P^ /v  4,  say.  For 
stimulated  Raman  scattering,  a  prominent  nonlinear  effect  in  waveguides. 
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one  obtains  z„  /z  This  value  represents  an  upper-bound 

1?  p 

estimated  for  z„/z  because  the  value  of  z  actually  depends  on  the 

t  p  P 

initial  density  of  spontaneous  photons;  moreover,  note  that  z^/Zp 

will  decrease  as  one  moves  to  longer  wavelengths  (see  Ref.  13).  Thus  the 

condition  z„^  z  will  be  satisfied  in  a  variety  of  circumstances, 

r  ^  p 

implying  that  self -focusing  should  be  observable  before  parametric 
effects.  Even  when  Zp^-^Zp  ,  parametric  processes  may  be 
distinguished  from  self -focusing  because  their  observation  requires 
overcoming  linear  losses  (i.e.,  g  P/A^#^,  where  is  the  waveguide  loss); 
self-focusing,  on  the  other  hand,  is  an  instability,  which  to  a  first 
approximation  is  insensitive  to  loss.  Various  other  parametric  processes 
may  be  suppresssed  or  distinguished  from  self -focusing  in  other  ways. 
Stimulated  Brillouin  scattering,  for  example,  can  be  suppressed  by 
controlling  the  dopant  profile  of  the  waveguide  and  in  extreme  cases 
could  be  eliminated  entirely  by  using  a  source  with  a  frequency  bandwidth 
larger  than  typical  acoustic  phonon  frequencies.  Multiple-wave 
parametric  processes  will  be  large  only  for  special  "phase-matched" 
frequencies;  moreover,  tight  tolerance  on  variations  in  diameter  and/or 
in  numerical  aperture  along  the  fiber  must  be  satisfied.  These 
considerations  serve  to  distinguish  self -focusing  from  other  nonlinear 
effects  operative  in  waveguides. 
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RADIATIOM  LOSS  IN  TAPERED  WAVEGUIDES 


1.  INTRODUCTION 

The  role  of  junctions  between  large-core  multimode  optical  waveguides 

and  small-core  single-mode  guides  in  integrated  optical  devices  has  led 

to  interest  in  tapering  between  fiber  waveguides.  This  interest  has 

12  3-7 

included  experimental  '  and  theoretical-numerical  studies  of 

tapered  structures.  An  important  theoretical  basis  for  several  of  these 
4-5 

studies  has  been  provided  by  the  step-approximation  method  developed 
by  Marcuse  '  for  tapers  in  step  waveguides.  This  approach  is  derived 
from  usual  perturbation  methods  by  regarding  arbitrary  waveguide 
deformations  as  a  succession  of  infinitely  many  infinitesimal  steps.  It 
has  provided  improvements  in  calculating  tapering  losses  beyond  the 
initial  development  of  the  perturbation  theory  of  small-wall 
distortions.^® 

We  use  an  alternative  approach,  the  eikonal  approximation 

method, for  the  calculation  of  the  total  radiation  loss  of  a 

guided  mode  that  is  due  to  a  symmetrically  narrowing  linear  taper.  This 

approach  has  been  developed  in  the  engineering  literature  for  optical 

dielectric  waveguides  as  the  concept  of  local  normal  modes  (LNM's)  or 
16-18  21 

tapered  modes,  '  paralleling  earlier  work  on  metallic 

waveguides.  However,  the  word  eikonal  has  also  been  used  by  physicists 
to  describe  a  similar  approximation  scheme  involving  a 

Wentzel-Kramers-Br illouin  (WKB)  approximation  for  the  initial  state  in  a 

13  19 

scattering  matrix  element.  ' 

We  present  the  eikonal  or  LNM  approach  by  beginning  with  first 
principles  for  TE  wave  propagation  in  optical  slab  waveguides. 
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indicating  both  the  small  parameter  of  the  method  and  its  region  of 
applicability.  The  radiation-loss  calculation  is  based  on  a 
weak-coupling  approximation  in  which  a  propagating  undepleted  symmetric 
pump  mode  couples  with  radiation  modes  in  the  tapering  region.  This 
represents  no  inherent  limitation  on  the  eikonal  (LNM)  method  but  does 
facilitate  comparison  with  perturbation  theory  and  with  the 
step-approximation  method. 

The  total  radiation-loss  result  is  reduced  to  similar  results  for  the 
step-approximation  method  by  keeping  only  the  more  important  of  two 
matrix  elements  on  which  the  coupling  theory  is  based.  Radiation-loss 
calculations  for  tapers  have  been  made  previously  by  using  other 
techniques, but  the  eikonal  (LNM)  method  has  not  been  directly 
e^.ployed  for  this  purpose. 

2.  RADIATION  LOSS  OF  A  SYMMETRICAL  LINEAR  TAPER 

From  the  time-independent  wave  equation  for  TE  modes  and  refractive 
index  n(r,z),  with  longitudinal  and  transverse  variables  z  and  r, 
respectively,  the  electric  field  is  expanded  in  terms  of  longitudinally 
dependent  coefficients  A^(z)  and  transverse  functions  (J>^(r,z) .  The 
latter  are  parametrically  dependent  on  z  and  are  defined  to  satisfy 
eigenvalue  equations  of  the  form 


41  Ct,*)  =  (*1 


(1) 


Where /2^(z)  is  the  corresponding  eigenvalue  that  is  also  assumed  to  be 
parametrically  dependent  on  z  and  =  w/c.  Use  of  the  expansion 

m. 
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in  the  wave  equation  and  mode  orthogonality  result  in  an  infinite  set  of 
coupled  equations 


^  $  *(a)  ^  2  F  A  =  0 

“  '^y\.  '  '  X  — '  *"*' - -  XVI.  v*v 

d  t  di 


(2) 


We  have  employed  the  shorthand  notation 


(2a) 


(2b) 


Equations  (2)  provided  the  basis  for  analysis  of  physical  processes  such 
as  reflection  amd  mode  coupling.  The  only  limitations  on  the  above 
theory  are  that  guide  changes  must  be  gradual  enough  so  that  the 
adiabatic  assumption  of  the  parametric  z  dependence  applies.  This  is 
equivalent  to  the  condition 

(3) 


where  >  is  the  wavelength  of  the  radiation  in  the  medium. 

In  order  to  apply  the  above  general  theory  to  our  tapering  problem, 
we  first  consider  the  propagation  of  an  initially  strong  mode  Ej^  = 
which  is  undepleted  on  interaction  with  the  other  modes.  Equation  (2) 
then  becomes 

+  f  +  *^<4  («l)  =  0 

Ai'  ^ 


which,  for 


(4b) 


and  has  the  solution 
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- 


'"(.I  t  " 


(4c) 


2  2  2 

where  terms  of  the  order  (d  k.^/dz  )  and  (dkj^/dz)  have  been 
neglected.  Equation  (4c)  is  valid  away  from  turning  points. 

A  reflection  coefficient  can  be  defined  for  pump  mode  Ej 
using  the  field  expansion  E  =Ej  +  E^  (2)  giving 


(a'irfa'  ,  ,  . 

-  +  (^)  ^  0  ^  I  '»)  1 

rf  T' 


(5) 


•a  ((9 

Here  k'^  (z)  =  d'fc^/di  ,  and  the  effect  of  the  neglected  second-order 

derivatives  for  A,  (z)  in  Eg.  (4c)  on  A,  (z)  has  been  retained  in  the 

right-hand  side  of  Eq.  (5).  In  Eq.  (5)  the  fast  and  slow  dependence  of 

A.  (z)  can  be  separated  out  by  focusing  on  a.  (z) : 

-tR  *K 


/■e 


Si  (e) 


(6) 


where 


^*1  =  j  («'l 


By  substituting  Eq.  (6)  into  Eq.  (5)  and  neglecting  the  second  derivative 
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of  (z)  in  a  slowly  varying  amplitude  approximation,  integration  of 

the  result  for  a  disturbance  located  between  0  and  z  gives  the  reflection 
coefficient 


2 

/  J \  e  di' 

R  - 

>> 

• 

1 

0 

A  similar  analysis  for  coupling  between  two  modes  with  dissimilar 
transverse  structure,  a  pump  mode,^  and  a  mode  n  with  forward-  and 
backward-moving  slowly  varying  components. 


results  in  forward  and  backward  and  coupling,  given  by 
^  Vi  V 

(a)  A-x 

The  mode  conversion  or  forward  scattering  of  a  propagating,  guided 
mode  to  radiation  modes  that  is  due  to  a  linear  taper  is  an  extension  of 
the  above.  The  geometry  of  this  problem  is  a  longitudinally  symmetric 
linear  narrowing  of  a  step-index  slab  guide  with  infinite  cladding.  The 
longitudinal  direction  is  divided  into  three  regions  for  which  the  guide 
half-width  L(z)  assumes  the  values  <L 
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= 


Ai* 


2.  <  ^ 


within  the  guide,  for  |j<l  ^  Ii(2),  the  index  of  refraction  is  n^^ 

whereas  outside,  for  |x  1  >  L(2),it  is  where  n^^  >  is 

assumed.  Note  that  AL  =  i —  -  L^  , 

For  TE  modes,  the  wave  equation  is  simply  (  /ix*  )+  (  ^ 

2  2 

+  n  k  E  =  0,  where  n  =  n,  n.,  in  the  appropriate  transverse 
o  i ,  <d, 

regions.  From  the  field  expansion  E(x,2)  =  X,  A  (  i  )  4*>n.(  )  , 

yw.  ' 

the  eigenvalue  equation  for  (x,2)  becomes 

t  k"  it/  (P^  (10) 

;  X* 

where  /3^(z)  is  equivalent  to  the  longitudinal  wave  number.  When  we 
define 


fxl  4  L 


[xl  >  L 


the  symmetric  normali2ed  (orthogonal)  eigenfunctions  take  the  form 


A/^  (*)  (el  K 


(12a) 


/\/  (e)  Ky^  (el  1^)  e. 


l<i  >L  (12b) 


~l  “  /z 

where  the  normalization  factor  is  N^(z)  =  [  L(z)  +  (z)l 

Since  the  linear  tapering  presents  a  symmetric  disturbance  to  the 

propagating  guided  modes,  it  provides  coupling  from  symmetric  modes  to 

other  symmetric  guided  modes  and  between  antisymmetric  and  other 

antisymmetric  guided  modes.  We  simplify  the  analysis  by  treating  only 

the  former  coupling.  Note  also  that  the  z  dependence  of  K^(z),  Y ■^iz) 

in  Eqs. (12)  is  due  to  the  matching  of  the  x  components  of  the  magnetic 

14 

field,  resulting  in  the  (symmetric  mode)  eigenvalue  equation 

£a)  £-(^)  =  ^wv  (*) 


If  we  consider  once  more  the  propagation  of  a  single  pump  mode  =A^ 
^^(x,2)j,  we  find  that  the  longitudinal  coefficients  of  this  mode 
satisfy  Eq.  (4a)  with  the  solution  from  Eq.  (4c): 


£*)  e. 


i  (»i 


(13) 


This  relation  is  valid  away  from  turning  points,  provided  that 

>4 £*)  '*'/  C. 


V 


It  will  prove  advantageous  to  discuss  the  coupling  of  the  above 
guided  mode  to  a  continuum  of  radiation  modes  by  using  discrete 
notation.  This  is  physically  equivalent  to  surrounding  the  tapered  step 
guide  with  an  ideal  conductor  at  x^^  t  R  assumed  to  be  large.  We  will 
recover  the  continuous  notation  as  the  latter  distance  becomes  infinite. 
The  TE  field  expansion  ^  K^(z)  ''P^(x,z)  in  the  wave 

equation,  together  with  the  wave-number  definitions 
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IX  I  L 


(14a) 


(>C|  >  L  (14b) 


leads  to  the  following  transverse  eigenfunction  structure  for  the 

14 

symmetric  radiation  modes  : 


(i) 


(8^  <r^  X 

•2.  [  [/‘f.  lx  I  "  (*)] 


IXI  4  L 

IXI  ^  L 


(15a) 


(15b) 


The  radiation  modes  are  assumed  to  exist  for  z  ^ Zq  -  A2  ,  and  the 
continuity  of  the  related  magnetic-field  components  at  the  transverse 
boundaries  gives 


>1^(8)  '  I  >1^  (J^) 


c 


Vi- 


Lfi,  r  ') 

Cto.  --  .^<r^  i.{il  I 

i  n  j 


(16) 


By  considering  the  slowly  varying  amplitudes  bq(z),  where  = 
b  exp  (  +  i/5,  a  ),  and  neglecting  the  second  derivative  of  b  (z)  ,  we 
can  obtain 


b  2.]  =  L  «=  -777 

'  »  -2  I  A 


1 

1 

r 

-J 

1 

1 

2. 

(S-  A  F7 

j 

xr  Djt 

1  1 

-  (fi' 

, 

-t/S.2' 

£  ^  c(2' 


(17) 


The  total  forward-scattered  radiation  loss  can  be  evaluated  from  Eq. 
(17)  in  the  limit  of  the  continuum  by 
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Total  radiation 


1 


It.,!' 


^  ‘  / ■k^_  ) 


mode  loss 

In  this  limit,  the  radiation  mode  wave  numbers  of  Eq.  (14)  become 


(18) 


A  (p\ 


'*■  Ir  ^ 

%  - 


A' 


(19a) 


—  *  0- 


If)  ^  (>’/  -  )  <^o*  +  =  'fc^ 


(19b) 


where  Eq.  (19b)  defines  K  .  Since  we  must  sum  over  all 

c 


forward -propagating  radiation  modes,  we  restrict  the  range  of  integration 
on  ^  to  0  ^  4  ”2*^0  avoid  imaginary  .  With  these,  Eq.  (18) 


becomes,  as  b 


b(^) 


Total  radiation  ^ ‘*\i '''•  |if(o)(^ 


mode  loss 


J 


Also,  it  can  be  shown  form  the  definition  [Eq.  (2a) 1  that 


1 


(19) 


-  /  ,  a  (i]  M.  re)  /  > 

~  -  (  K  ^  O'*  L  Ci] 

U\  -  fr.  ^  O?  I  V 


(20) 


Furthermore,  F^j^Cz)  of  E<3*  (2b)  provides  higher -order  corrections 

in  Af:  to  G.^(z)  and  hence  may  be  neglected  in  the  following 
Ai  ' 

first-order  calculation.  Further  approximating  ■)^^(z)  by  /S'^(z)  from  Eq. 
(4b) ,  we  obtain  our  main  result  after  some  algebra 


Total  fractional  forward -scattered 


radiation  loss 


(m) 


i 


a 


e-'l  /  t  <  .  V..*rLV'* 


(21) 
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For  convenience,  we  summarize  the  locations  of  the  definitions  of  the 
symbols  of  the  fractional  loss,  Eq.  (21);  are  defined  by 

Eqs.  (ll);p(/»),  <r  {f>) ,  ^  are  related  by  Eqs.  (19);  K^=  (n^-nj) 
is  defined  in  Eq.  (19b);  L(z),  A  z,  are  defined  in  Eq.  (9). 

Also  (z)  is  found  after  obtaining  (z)  from  the  eigenvalue  equation 
with  L  =L(z)  from  Eq.  (9).  It  is  only  slightly  different  from^^  ,  which 
is  obtained  by  using  L=L-. 

The  fractional  forward-scattering  radiation -loss  result  of  Eq.  (21) 
has  been  limited  only  by  the  first-order  approximation,  including  only 
However,  the  relatively  mild  approximations  involved  are  still 
insufficient  for  analytical  evaluation.  We  may  obtain  the  main  effects 
by  approximating  the  integrand  of  Eq.  (21)  to  zero  order  in  the 
longitudinal  dependence  for  the  amplitude  and  to  first  order  in  the 
phase.  This  would  be  reasonably  accurate  for  weak  tapers.  The  essence 
of  the  eikonal  approach  is  that  it  correctly  takes  into  account  overall 
phase  changes  of  the  wave  while  propagating  through  a  disturbance.  This 
is  the  main  distinction  of  the  eikonal  approach  as  compared  with 
perturbation  theory. 

Furthermore,  we  expect  the  strongest  coupling  to  the  radiation  modes 
to  come  from  higher-order  guided  modes  since  the  latter  are  less  confined 
to  the  core  and  have  stronger  cladding  fields.  For  such  guided  modes, 
the  value  of  K^(z)  L(z)  can  be  seen  to  lie  near  Ztt  from  a  graphical 
evaluation  of  the  eigenvalue  equation.  With  this,  the  zero-order  value 
for  K^(z)  is  l'(z) /ttw and  the  values  are  similar  for  the  other 
variables  related  to  (z)  from  Eqs.  (11).  The  lower-order  guided 


modes,  which  are  more  tightly  confined,  have  K^(z)  near  (21+1) /2 

I  TT  /  L  (a)  )  !s  (21+1) /2  (  '*’'/l_),  and  hence  the  matrix  element  G,  Az)  of 

r* 

Eq.  (20  )  becomes  small.  For  higher-order  guided  modes,  to  zero  order 
in  the  z  variation. 


0 


Xtr  I 


1__ 


(22a) 


(22b) 


(22c) 


Recall  that  K^_  is  the  value  of  Kj((2)  from  the  eigenvalue  equation 
evaluated  for  L(z)  =  L-.  Also,  in  the  denominator  of  the  integrand  in 
Eq.  (21),  the  following  rather  complicated  factor  will  be  approximated  as 


Ctm.  (rt~ 


'u  r  ,  r  ~\  'l>- 

=  o- oL\~  tv  ~  It  !  (22d) 
L  L.  2.  i 


where  we  have  neglected  the  oscillation  that  is  due  to  cos^  0~  L(z)  over 
the  transition  region. 

We  now  use  the  guide  width  variation  to  define  the  variable 


it?\  ^  Ly.  li) 

At 


when  the  above  approximations  are  combined,  Eq.  (21)  becomes 
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Total  fractional  ''?■** 
forward-scattered  /  ^ 


•u 


radiation  loss  = 


/r^i(  dp 


(23) 


where  -Cipf <■  o-L^  i' +  <■  f  /Sj^  1 1-  ^ 

L  ^'4*  fix,*- 


i'M  ^ 


^  Cr<-_ 

+  €■ 


-i(rL_o(  (i'l-ifii.'-*-  M  /^/o  (  I  -  o<(8"l] 

i -M  /Sj^ 


The  resulting  quadratic  z  dependence  in  the  exponentials  in  Eq(23) 
permits  the  z  integration  to  be  expressed  in  terms  of  Fresnel  integrals 


r  X  r  K 

5"  (x  ]  =  /  JITw^^^W  =  l  c-tna,  JL  dw 

\  ^  ‘  j  ^ 

O 


(24) 


With  these,  the  result  of  the  z  integration  in  Eq.  (23)  is 

'  fJ:  [■^''’'1  '  I  [AcS  *- *  lAt^'-^(AK]‘ 

•f  ^  ACi^^Cy  "i  ^  ^  /Q-t'vu  )  ^  ^  Cy  — 


(25) 


Here 


AC,  =  C(u.j  r.  fvJ^Cfv,>| 


A  S,  =  -  Ifu/  A  S-y  =  5  -  r  (v  I 


(25a) 


and  we  have  defined  the  following  dimensionless  parameters  as 


I  /J,  A2  At 

I  <0  — 


(25b) 
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f  (^u  ~41p] 
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(25c) 


(25d) 


By  using  Eqs.  (25)  in  Eq.(23),  we  obtain 
Total  forward-scattered 


(26) 

where  |  -f-(^lj  denotes  the  ^  dependence  within  the  Fresnel  integrals  as 
indicated  by  the  part  of  Eq.(25  )  within  the  braces.  Because  of  the 
nature  of  the  approximation  in  Eq.  (22a)  ,Eq.  (26)  as  written  is  limited 
to  higher-order  pump  modes.  However,  if  the  dependence  on 
is  changed  to  K.  Y.  ,in  Eqs.  (25)  and  (26),  the  latter  are  also 
suitable  for  describing  the  radiation  loss  of  the  lower-order  modes. 

The  graph  in  Fig.  1  illustrates  the  fractional  radiation  loss  for 
four  cases  of  taper  step  size  using  Eq.  (26)  forA^=l>i~\  free-space  guided 
mode  wavelength.  The  four  cases,  Au  »  7. 5, 5, 2. 5,  and  1  urn,  correspond  to 
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0  0  0  ^ 

tapering  angles  of  approximately  8.5,  6,  3,  and  1  ,  respectively,  for  an 
initial  guide  half-width  of  L-  =25  and  a  tapering  region  length  Az  =50 
/um.  The  refractive-index  values  chosen  were  nj^=1.05  and  n2=l-00. 

These  parameter  values  permit  only  il  of  the  15  modes  that  propagate  in 
the  initial  wider  guide  to  propagate  through  all  four  tapers. 

These  four  cases  show  a  rapid  increase  in  total  radiative  loss  as  the 
mode  number  approaches  the  cuttoff  value,  par ticurlarly  for  the  smaller 
taper  angles.  This  tendency  is  in  agreement  with  similar  results 
observed  in  Ref. 5.  Note  that  the  near-linear  behavior  of  the  total 
radiation  loss  on  a  logarithmic  scale  for  the  smaller  angles  indicates  a 

(caturji 

roughly  exponential  dependence  on  mode  number  jL  ,  RAD~«-  .  Analogous 
results  are  seen  in  Fig.  2  in  which  the  same  parameters  are  used  for 
plots  of  radiation  versus  step  size  for  the  modes  indicated.  Again 
RAD ~e  .  These  trends  could  serve  as  useful  rules  in  making 
comparative  estimates  of  the  radiative  loss  in  such  structures.  Of 
course  the  accuracy  of  this  approach  decreases  as  both  ai-  and  the  mode 
number^  increase  since  they  both  contribute  to  an  increase  in  the 
left-hand  side  of  Eq.  (3).  Similar  plots  obtained  using  perturbation 
theory  indicated  higher  total  radiation  loss  for  the  higher  modes  as 
compared  with  the  plots  of  Fig.  1.  Modes  beyond  those  shown  are  totally 
reflected  and  must  be  treated  by  a  different  analysis  because  of  the 
breakdown  of  the  WKB  approximation  inherent  in  the  theory  through  Eqs.(4). 
3.  DISCUSSION 

The  suitability  of  the  eikonal  or  LNM  prodecure  to  tapering  problems 
has  suggested  the  main  calculation  of  this  section,  a  derivation  of  the 
total  radiation  loss  of  a  guided  mode  that  is  due  to  a  symmetrical 
linear  taper.  We  have  limited  our  results  to  a  weak -mode -coupling 
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theory  in  which  an  undepleted  pump  mode  transfers  energy  to  radiation 
modes  and  the  coupling  equation^ are  not  solved  in  a  self-consistent 
manner. 

An  evaluation  of  Eq.  (21)  with  zero-order  amplitude  and  first-order 
phase  approximations  leads  to  Eq.  (26),  which  involves  Fresnel 
integrals.  This  shows  the  total  radiation  loss  to  be  essentially 
proportional  to  the  gradient  (  Ai-  /  ) .  Analytical  evaluation  of  the 

final  result  of  the  integrations  of  Eq. (21)  is  difficult,  but  it  would 
serve  as  an  appropriate  starting  point  for  further  numerical  studies. 

In  contrast,  the  approximate  expression  to  Eq.  (21),  Eq.  (26),  requires 
only  a  relatively  simple  one -dimensional  evaluation.  The  plots  obtained 
numerically  from  Eq.  (26)  illustrate  the  application  of  the  theory  and 
provide  quantitative  results  for  the  total  fractional  radiative  coupling 
of  guided  modes  in  the  undepleted  pump  approximation.  As  expected.  Pigs. 
1  and  2  show  a  strong  increase  in  radiative  loss  for  increasing  At.  and 
for  modes  approaching  cutoff.  The  approximate  exponential  dependences 
noted  in  the  figures  have  not  been  discussed  before.  Furthermore,  the 
total  radiation  loss  for  tapers  has  not  been  previously  calculated  by 
using  the  eikonal  or  LNM  method,  and  this  quantity  particularly 
facilitates  comparison  with  experiment. 

The  simplicity  and  general  applicability  of  the  eikonal  or  LNM  method 
under  the  proviso  of  small  gradients  have  been  indicated  in  the 
development  of  the  theory.  The  approach  is  particularly  well  suited  to 
investigating  tapered  structures  but  is  also  readily  adaptable  to  any 
gradual  longitudinal  disturbance. 
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PIGU1<.;  CAPTIONS 


FIG.  I  Total  fractional  radiation  loss  vs.  mode  number  from 
eq.  (^(p)  for  radiation  •»  J  ^  in  a  slab  waveguide  of  initial 
halfwidth  L_  =  25^4  and  tapering  region  Az  =  50^  , tapering  angle 
e  =  —  ,  n,  =  1.05,  n^  =  1.00:  (a)  &  i~  =  7.5^  (b)  A  L  =  5^ 

{c)  At-  =  2.5^  (d)  At  =  1^  .  An  approximate  exponential  dependence 
on  mode  number  is  noted  for  Al.  =  1  . 

FIG.  2  Total  fractional  radiation  loss  vs  taper  step  size 
from  eq.  (a-fc)  for  the  same  parameters  as  in  Fig.  1.  (a)  mode  ^ 

(b)  mode  H  (c)  mode  1  (d)  mode  o. 


-52- 


K I  cunt;  1 


LOG,^  (TOTAL  FRACTIONAL  RADIATION  LOSS) 
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